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ABSTRACT 

Aims. We present a method that allows us to estimate the distance from the continuum source located in the center of AGN to the 
highly ionized gas called warm absorber. 

Methods. We compute a set of constant total pressure photoionization models compatible with the warm absorber conditions, where 
a metal-rich gas is irradiated by a continuum in the form of a double power-law. The first power-law is hard up to 100 keV and repre- 
sents a radiation from an X-ray source, while the second power-law extends downwards from several eVs and illustrates a radiation 
from an accretion disk. 

Results. When the ionized continuum is dominated by the soft component, the warm absorber is heated by free-free absorption, in- 
stead of Comptonization, and the transmitted spectra show different absorption line characteristics for different values of the hydrogen 
number density at the cloud illuminated surface. 

Conclusions. This fact results in the possibility of deriving the number density at the cloud illuminated side from observations and 
hence the distance to the warm absorber. 
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1. Introduction 

Many Active Galactic Nuclei (AGN) exhibit numerous absorp- 
tion features of highly ionized material in their UV/X-ray spec- 
tra. Such gas, called Warm Absorber (hereafter WA), is located 
on the line-of-sight towards the observer and it is illuminated by 
radiation originating from the active nucleus. 

Most of the observed absorption lines are blushifted, sug- 
gesting that the metal-rich gas is outflowing. We still don't know 
how such a wind is powered and where it is lunched, neverthe- 
less, due to high resolution spectroscopic observations of AGN 
from FUSE, Chandra, XMM and other satellites, we are able to 
make some d iagnostics on the physical condi t ions of the WA ( for 
a review see: lKrissll2004t iBlustin et al.ll200H iRozanska 2007]). 

The observed absorption lines velo city shifts are of th e order 
of 10 4 km s" 1 in case of UV absorbers dGabel et alj2003l) and of 
the order of 10 2 to 10 3 km s _1 in the case of X-rays ( Kaspi et"ai1 
l200lHKaastra et al.l2002h . The column density of the WA is gen- 
erally estimated to be about 10 21 ~ 23 ctrT 2 . The absorbing gas 
comprises different ionization phases, corresponding to temper- 
atures from about 10 5 K, when iron is partially ion ized, up to a 
few 10 7 K, when iron is almost completely ionized |Netzer et al.l 
120031: IKrongold et al.|[2003t ISteenbrugge et al.ll2005ir From the 
above hints, photoionization modeling of the WA is made, trying 
to answer to the most important and unsolved question: how far 
from the continuum source is the absorbing gas located? 

However, there is one difficulty which prevents us from an- 
swering this question. It is well known that photoionized models 
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of a cloud illuminated by a single X-ray power-law, typically 
ranging from 0.01 keV up to 100 keV (her eafter we will cal l this 
a hard X-ray illumination) are degenerate dRozaf iska 2007). We 
cannot distinguish between clouds with the same ionization pa- 
rameter, but different hydrogen number density no, at the cloud 
illuminated surface, and at different locations. The transmission 
spectrum from a rarefied cloud (no ~ 10 5 cm 4 ) located farther 
away i.e. at ~ 0.1 pc, is identical to the spectrum of a dense 
cloud («o ~ 10 10 ctrT 3 ), situated at ~ 0.0001 pc from a contin- 
uum source. The first case is consistent with the WA being co- 
spatial with the NLR (narrow line region) or dusty torus, while 
the second case corresponds to the closest neighbourhood of an 
accretion disk and BLR (broad line region). 

Some estimations on the distance to the WA were done us - 
ing variability studies dNetzer et al.ll2003l: IKrongold et al.l2005l) . 
However, us i ng the same 900 ksec Chandra data of NGC 3783, 
iNetzer et al.l (|2003) found upper limits for the location of three 
ionization phases to be 3.2, 0.63, and 0,18 pc respectively, while 
IKrongold et al.l d2005l) estimated the location of the high ioniza- 
tion phase at 0.0029 pc and low ionization phase at 0.0004-0008 
pc. 

In this paper we show that the degeneracy of photoionized 
models breaks down when the WA is illuminated by an AGN 
continuum including the thermal disk radiation, represented here 
as a second power-law component. The disk contribution is par- 
ticularly important in the case of quasars, since their broad band 
spectra are clearly dominated by the soft disk emission. In this 
article we explain why double power-law models are not degen- 
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erate and claim that using such models provides a way to deter- 
mine the distance to the WA in some objects. 

The structure of the paper is as follows: In Sec. [2] we ex- 
plain how to estimate the distance to the WA from photoionized 
models. Section [3] describes the photoionization models used in 
this work. Results for different spectral shapes are presented in 
Sec.|4] Discussion of the models is presented in Sec. [5] and main 
conclusions of our work are drown in Section [6] 

2. The Distance to the Warm Absorber 

The distance from the central engine to the WA is a clue issue 
for understanding the wind geometry. We can derive this dis- 
tance from observations using photoionized models, which are 
parametrized by the so called ionization parameter, defined at the 
cloud surface. There are three ionization parameters being used 
by different groups and photoionization codes, but the idea is the 
same: to define the amount of ionized radiation which reaches 
the cloud. Here, we discuss the problem of distance derivation 
using the ionization parameter defined as: 



where L I0 „ is the ionizing source luminosity, n is the hydro- 
gen number density, and R is the distance from an illuminating 
source. Nevertheless, the problem of distance derivation does not 
depend on the ionization parameter used by different groups, and 
refers to the value of this parameter on the cloud surface i.e. for 
intrinsic luminosity, and for surface hydrogen number density, 
hereafter «o- 

From observations, we are able to get the luminosity of the 
object as a separate measurement, and the ionization parameter 
on the cloud surface is obtained by fitting the data with a given 
photoionization model. But we are still left with the noR 2 prod- 
uct, meaning the distance to the WA depends on the value of the 
number density at the cloud illuminated surface. 

Below, we show how to determine the number density at 
the illuminated side of the WA using photoionization modeling. 
Then, after simple inversion of Eq. [2] we obtain the distance to 
the absorber. 

3. Photoionization modeling 

In this paper all photoioniz ation models a re computed using 
the tit an code developed by iDumont et al.l d2000t) ; ICollin et"aH 
titan is well-suited both for the study of optically thick 
and thin media, such as the WA. It computes the gas structure 
in thermal and ionization equilibrium assuming non-LTE, and it 
provides the reflected, emitted outward, and absorbed spectra. 
titan treats the transfer of both the lines and the continuum us- 
ing the ALI (Accelerated Lambda Iteration) method, which pre- 
cisely computes line and continuum fluxes in a self-consistent 
way. Our atomic data include ~ 10 3 lines from ions and atoms 
of H, He, C, N, O, Ne, Mg, Si, S, and Fe, all elements having 
cosmic abundances. 

The total (gas+radiation) pressure is computed after the 
plasma temperature has been determined by the thermal balance 
equation. The pressure is kept constant within the medium thus 
allowing to determine the density profile. The model is iterated 
until convergence. We have shown that the assumption of con- 
stant pressure al lows the ionized gas to be naturally stratified due 
to illumination (Ro zanska et al.l2006HGoncalves et al.l2006l) ; as 
a result, the cloud comprises different ionization stages which is 
consistent with observations. 



In all our models we assume a plane-parallel geometry, with 
the WA medium being illuminated on one side by a given in- 
cident continuum. We assume the external radiation to hit the 
cloud perpendicularly to its surface. 

In the frame of photoionization modeling of the warm ab- 
sorber, there is no possibili ty to reconstruct velocity structure ob- 
served in some cases (e.g. ISteenbrugge et al.ll2005l) .where each 
ion has slightly different velocity. We can only consider the case 
where all ions have the same turbulent velocity, and it is taken 
into account in our code. 

The models are parametrized by the ionization parameter § 
at the cloud surface, the total column density Nh, the hydrogen 
number density at the cloud illuminated surface «o, and the inci- 
dent flux Spectral Energy Distribution (SED). 

The parameter i; is very well constrained on the surface of an 
illuminated medium. But if we want to trace ionization proper- 
ties deeper in the gas, we use the so called dynamical ionization 
parameter defined as: 

rn ^ Lj on 1 E ion Prod 

a ~ AnckT ~ 4ncR 2 ~riT ~ cP gas ~ ~P^~ S ' 

where c is a velocity of light, T is temperature of the gas, F, o; , 
is a flux which affects the cloud, and P ra d and P gas are radiation 
and gas pressure respectively. Below, we describe the ionization 
structure of presented clouds using relation E(T). 

In this paper we consider two possible shapes of SED first: 
a single power-law roughly representing the spectra of Seyfert 
galaxies, and second: a double power-law representing the spec- 
tra of disk dominated bright quasars. In both cases, the power- 
laws are modified by exponential high- and low-energy cut-offs 
according to the formula: 

F E = E' a exp(-E/E max )exp(-E mi „/E), (3) 

where F is the flux, E - energy, a - spectral index, £„„„ and 
E max are low and high spectral ranges respectively. In the case of 
double power-law models, we have additional parameter, which 
is the relative normalization of both power-law components, i.e. 
the ratio of their integrated fluxes denoted as FxjF so f t . This rel- 
ative normalization and exponential cut-offs ensure that there is 
no discontinuity at the point, where two power-laws merge. 

4. Results for the different shape of illuminated 
continuum 

To illustrate the determination of distance by photoionized mod- 
els, we consider three clouds of different hydrogen density num- 
ber at the illuminated surface: «o = 10 6 , 10 s , and 10 10 cm -3 . 
The ionization parameter on the cloud surface is the same in the 
three cases and | = 10 5 . The total column density of each cloud 
is taken to be close to the maximum column density for which 
the temperature profile is thermally stab le for a given spectral 
energy distribution (see Rozanska 2007|). In all cases, the total 
column density is of the order of 10 23 cm -2 , and we assume the 
turbulent velocity of the gas equaled 300 km/s. 

For a single power-law irradiation, the primary continuum 
spectrum ranges from 0.01 up to 100 keV and has a spectral 
index equal to 0.8. In the case of a double power-law, the softer 
power-law ranges from 10~ 7 up to 0.0136 keV with a spectral in- 
dex a so f t = 0.3, and the hard continuum component has a spec- 
tral index ax = 0.8, ranging from 0.0136 up to 100 keV. We 
consider three cases of relative normalization of both power-law 
components : (i) Fx/F so f t = 0.01, corresponding to a domina- 
tion by the soft component; (ii) Fx/F so f t = 0.5, typical of an 
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Fig. 1. Ionization structure (a), temperature structure (b), and transmitted spectrum (c) for a WA illuminated by a hard single 
power-law continuum for three different values of Hq: 10 6 crrT 3 (green line), 10 8 ctrT 3 (red line), and 10 10 cirT 3 (blue line). 




Fig. 2. Ionization structures (a), temperature structures (b), and transmitted spectra (c) for the WA illuminated by double power-law 
for three different values of «o : 10 6 cm -3 (green line), 10 8 cm -3 (red line), and 10 10 crrT 3 (blue line). Upper panels correspond to 
Fx/F s „f, equal 0.01, middle panels to 0.5, and bottom panels to 8. 



intermediate case; and (iii) Fx/F so f t — 8, where the hard com- 
ponent dominates. 

4.1. Hard, single power-law photoionized models 

The comparison of the clouds with different no values irradiated 
by a hard, single power-law is presented in Fig.Q] Panel (a) rep- 
resents the ionization structure through the WA, which is iden- 
tical for all no values. The maximum temperature on the illumi- 
nated side of the clouds is quite high, of the order of 4 x 10 7 K. 
Going deeper into the WA, the temperature structure displays 
a strong gradient (panel (b)), which is consistent with the gas 
structure bei ng in hydrostatic equilibrium and illuminated by 
hard X-rays dRozariska et al.ll2002tlMadej & R6zafiskall2004l) . 

The transmitted spectra are presented in Fig. Q] (c). They 
differ only in normalization, which means that the low density 
cloud, with no = 10 6 cirT 3 , located further away from the con- 



tinuum source (see Eq. 0, has the same spectral features as the 
dense cloud of no = 10 cm 4 , located two orders of magnitude 
closer to the source. The detailed comparison of lines for those 
two cases is presented in Fig. |5](bottom panels). Lines are iden- 
tical in all spectral bands, therefore, models with a single, hard 
power-law illumination are highly degenerate and useless for the 
purpose of determining the distance to the WA. 



4.2. Double power-law photoionized models 

The situation changes when we consider a primary continuum 
described by a double power-law dominated by its soft compo- 
nent. In Fig.|2]we show results for clouds illuminated by a dou- 
ble power-law with increasing Fx/F so f t ratio, from 0.01 (upper 
panel) up to 8 (bottom panel). When the continuum soft compo- 
nent dominates, the ionization and temperature structures of the 
considered clouds differ for different values of «o- Separation be- 



4 



A. Rozariska et al.: Distance to the Warm Absorber in AGN 



tween models is clearer for lower values of the Fx/F so f t ratio. 
The equilibrium temperature at the illuminated face of the cloud 
is about one order of magnitude lower than in the case of a single 
power law illumination. 

For the highest value of no = 10 10 cm" 3 , the continuum spec- 
trum displays a significant absorption in the low energy band due 
to free-free absorption. Looking closer at the transmitted spectra, 
there is a big difference in the equivalent width (EW) of absorp- 
tion lines between clouds of two extreme values of no, and for 
FxlF so f t = 0.01; this is illustrated in Fig. |5]upper panels. For 
instance, EWs of HOg and H0 f lines for the case of «o = 10 10 
cm" 3 are ten times smaller than for the case of no = 10 6 cm" 3 . 
In X-ray domain, EWs of exampled lines: Fel9i, Fe20^, Fel7/; 
and Fel9// are respectively seven, five, fourteen and four times 
larger for the cloud with no = 10 10 cm" 3 than for the cloud with 
«o = 10 6 cm" 3 . Therefore, from the ratios of EWs of absorp- 
tion lines alone, we can select one particular photoionized model 
which fits the observations, and then having obtained the num- 
ber density at the illuminated cloud surface, we can derive the 
distance to the WA. 

The observed differences of photoionized models can be 
m uch higher for dif ferent spectral energy distributions as shown 
bv lRozanskal d2007t) . Below, we describe what causes this differ- 
ences in the case of a WA. 



5. Discussion 

In order to understand why single power-law models are degen- 
erate, and why after including a second illuminating component 
models can be effectively separated, we conducted an analysis 
of the basic physical conditions inside the clouds. We followed 
the behaviour of the parameters responsible for the hydrostatic 
and the ionized thermal equilibrium of the gas, i.e. pressure and 
efficiency of radiative processes. 
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Fig. 3. Comparison of the gas pressure (solid line), the radiation 
pressure (dashed line), and the total, i.e. gas+radiation, pressure 
(dotted line), for all computed models. The upper panel cor- 
responds to clouds with, «o = 10 6 cm" 3 , the middle panel to 
«o = 10 s cm" 3 , and the bottom panel to «o = 10 10 cm -3 . 



5.1. Pressure balance through the cloud 

All computed cloud models are dominated by radiation pressure, 
as shown in Fig. [3] For all cases, the radiation pressure (dashed 
lines) practically equals the total (gas + radiation) pressure (dot- 
ted line). There is nothing unusual in the behaviour of the gas 
pressure (solid lines) for all computed models. For a given «q 
the gas pressure is highest for the highest Fx/F so f t . Obviously, 
the value of the gas pressure on the cloud surface is proportional 
to the value of the hydrogen number density on the cloud illumi- 
nated side. 



5.2. Energy balance within the cloud 

The titan code solves the energy balance by calculating ra- 
diation gains and losses at each point of the gas slab. Fig. [4] 
shows the net heating minus cooling rate integrated over the 
whole cloud versus the ratio of the two continuum components, 
Fx/Fsoft. Each important mechanism was analyzed separately 
to identify its contribution to the total energy balance in the WA. 
Net Comptonization is represented in green, ionization minus 
recombination in blue, bound-bound heating minus cooling in 
black, and finally, free-free balance in red. 

The heating and cooling rates cover a significant range, so we 
decided to plot the results on a logarithmic scale. We arbitrarily 
assign a minus sign to those cases where total gains are negative, 
i.e. when cooling dominates (bottom panel). The sign is positive 
when heating dominates (upper panel). There is an arbitrarily 



discontinuity between the two panels in Fig. @] which does not 
have any physical meaning. A different type of line is used to de- 
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line), no = 10 cm" (dotted line), and no = 10 cm" (dashed 
line). Magenta hexagons represent the net heating minus cooling 
for the case of a single power-law illumination, they are obvious 
limits of the radiation rates when increasing the Fx/F so ft ratio. 

Hard X-ray radiation always heats the WA, even if the X-ray 
component is very weak. The same is true for the line cooling, 
which most probably is balanced by photoionization. The sit- 
uation changes dramatically for Comptonization and free-free 
processes. For Fx/F so ft decreasing from 1 and 0.1, Compton 
heating decreases and the WA gas starts to be even Compton 
cooled by soft disk photons. The same happens in the case 
of bremsstrahlung free-free processes. Bremsstrahlung heating 
starts to dominate when Compton heating becomes unimpor- 
tant. This means that, in the case of a WA illuminated by a dou- 
ble power-law, the presence of a soft, disk component changes 
the mechanism responsible for the high temperature equilibrium 
on the irradiated side of the cloud. Therefore, in such a situa- 
tion, clouds are heated mostly by bremsstrahlung, and not by 
Comptonization, like it is in case of clouds illuminated by a sin- 
gle, hard X-ray continuum. 
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Fig. 5. Comparison of spectral absorption lines for three energy bands: UV (left panels), Soft X-ray (middle panels), and Hard X-ray 
lines (right panels). Upper panels show the case where the WA is illuminated by a double power-law with Fx/F so f t = 0.01, while 
bottom panels correspond to irradiation by a hard, single power-law. Blue lines depict dense clouds of no = 10 10 crrT 3 , and green 
lines less dense clouds of no = 10 6 cm 4 . 



5.3. Absorption lines from double power-law models 

In Fig. |5j we focus on the absorption lines present in three 
spectral bands for clouds with no = 10 10 crrT 3 (blue line) and 
no = 10 6 crrT 3 (green line). The upper panels represent models 
with Fx/F so f t = 0.01, while the bottom panels depict models of 
clouds illuminated by a single, hard X-ray power-law. 

For strongly degenerate models, the absorption line spectra 
are identical and there is no possibility to distinguish them based 
on the number density «o just by fitting photoionization models 
to the observations. The situation is very different in the case of 
a double power-law illumination. The equivalent widths of some 
absorption lines differ for different no, giving the possibility of 
fitting one single photoionization model to the observed spec- 
trum. Therefore, having a grid of photoionized models providing 
the main absorption line ratios, we can estimate the number den- 
sity at the cloud illuminated surface, and thus the distance to the 
WA. Indeed, even modeling only the X-ray data of any AGN, we 
have to take into account that the WA is affected not only by an 
X-ray power-law from the corona, but also by a soft component 
from the disk. This is accounted for in this article and justifies 
the use of a double power-law extending from 10~ 7 to 100 keV. 



6. Conclusions 

Photoionized models of a WA usually assume that the illumi- 
nating continuum has the shape of a single power-law with the 
spectral index derived from X-ray observations of a particular 



AGN. In this paper, we have shown that it is important to include 
the broad-band continuum from the active nucleus to achieve a 
proper modeling of the transmission spectrum through the ion- 
ized WA gas. 

We have computed a set of models assuming illumination of 
the metal rich material by a double power-law incident contin- 
uum, where one component mimics the emission from an accre- 
tion disk, and the second component represents the hard X-ray 
emission coming most probably from a hot corona, or magnetic 
flare. This allows to break the degeneracy generally observed in 
photoionization models using a single hard X-ray illumination. 
By breaking this degeneracy we observe that models display 
different ionization and temperature structures for clouds with 
different number densities at the illuminated surface. Therefore, 
transmitted spectra look different and in principle it should be 
possible to identify the photoionized model which better fits a 
given observation. 

We have shown that this degeneracy breaks down due to 
a switch in the radiative mechanism heating the ionized cloud 
when Fx/F so f, decreases. For a single power-law irradiation, 
hard photons mostly participate to the Compton heating, and this 
process becomes responsible for the hot equilibrium temperature 
of the illuminated layers. When Fx/F so f, decreases, i.e. when 
the soft component starts to dominate, Compton heating is not 
efficient anymore, since the cloud begins to be more efficiently 
cooled by Comptonization, and a hot temperature equilibrium is 
established due to bremsstrahlung absorption. 
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Fig. 4. Net, i.e. Heating - Cooling gains of energy integrated over 
the cloud versus the ratio of hard X-rays to the soft power-law 
component. Net gains for Comptonization are represented by 
green lines, bound-free processes, i.e. ionization-recombination, 
are represented by blue lines, bound-bound processes, i.e. line 
heating-cooling, by black lines, and free-free processes by red 
lines. Axis are logarithmic, therefore we arbitrarily assigned a 
minus sign to those cases where total gains are negative, i.e. 
cooling dominates over heating (bottom panel). There is a dis- 
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ing (upper panel). For each process, solid lines represents clouds 
with «o = 10 10 cm 4 , dotted lines «o = 10 8 cirT 3 , and dashed 
lines «o = 10 6 cm 4 . Magenta hexagons mark gains for the case 
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In conclusion, by taking into account broad band illumi- 
nation from the active nucleus it is possible to determine the 
number density at the illuminated side of WA and thus its dis- 
tance from the irradiating continuum source. This method is 
complementary to t he m ethod based on variabi lity presented by 
iNetzer et all d2003l) and iKrongold et al.l d2005h . and can be ap- 
plied to bright quasars with a strong disk component. 
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